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Abstract 

Consider a pyramid made out of unit cubes arranged in square horizontal 
layers, with a ledge of one cube's length around the perimeter of each layer. 
For any natural number k, we can count the number of ways of choosing k 
unit cubes from the pyramid such that no two cubes are in the same horizon- 
tal layer; we can also count the number of ways of choosing k unit cubes 
from the pyramid such that no two cubes come from the same vertical slice 
(taken parallel to a fixed edge of the pyramid) or from two adjacent slices. 
Djakov and Mityagin first established, using functional analysis, that these 
two quantities are always equal (the enumerative interpretation given here 
is due to Thomas Kalinowski). Don Zagier supplied the first combinatorial 
proof of this result. We provide a new, more natural combinatorial proof. 

The authors of [O conjectured ([fl] Conjecture 1]; []2] |3] Theorem 4]) that the 
following identity holds for all positive integers n and integers k: 

E U? = E W^+^-j) a) 

where [n] := {1, 2, 3, . . . , n) as usual, (^') denotes the collection of /c-subsets 
of [n] all of whose elements are congruent to n modulo 2, and denotes 
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the collection of fc-subsets of [n] which do not contain two consecutive integers. 
Don Zagier gave a combinatorial proof of this result in the appendix to [[2l|3]|. It 
turned out that a few years prior to this, the same identity had been established, 
with a less direct, non-combinatorial proof, in the context of functional analysis 
(IS 13). In the present note, we establish the identity in the most straightforward 
fashion yet, by calculating the characteristic polynomial of a certain matrix in two 
different ways. 
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where the entries in the diagonal below the main diagonal are ai, a 2 , . . . , a n , the 
entries in the diagonal above the main diagonal are b\, hi, ■ ■ ■ , b n , and the remain- 
ing entries are 0. Let Xn{%) = det(xl — M) be the characteristic polynomial of 
M n in the variable x. We claim that for all n > 1 and all k G Z, the coefficient 

d kiU of x n+1 ~ 2k in Xn( x ) satisfies 
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One way to see this is by observing that both sides of Equation [3] satisfy the recur- 
rence 

r k,n = r k,n-l — a nb n r fc-l,n-2 (4) 

for n > 3, and they agree for n = 1,2. For the left-hand side of Equation |3] 
the recursion follows by expanding the determinant about the last row. For the 
right-hand side, it follows by considering the terms of total degree k with respect 
to the variables xi, x%, . . . , x n in the product 



(5) 
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after eliminating terms in which two consecutive indices appear. 

Specializing to the case a, = % and hi = n + 1 — i, the matrix M n becomes 
what Taussky and Todd term the "Kac matrix" S n (see iH for the long and im- 
pressive history of this matrixQ; though named for Victor Kac, the notation honors 
Sylvester, who found its characteristic polynomial in the 1850s). We only need 
the eigenvalues of S n ; these are computed, e.g., in [|6), but we present an original 
derivation of the spectrum in what follows. 

A non-zero vector v = [i>i i>2 • • • v n+ i] is an eigenvector of S n with eigen- 
value A if and only if the components Vi of v satisfy the relations 

ivi + (n - i)v i+2 = Xv i+1 (6) 

for all i £ {0, 1,2, ... , n}; notice that the coefficients of v and v n+2 will be 0, so 
this really does work. 

Let d be an integer between and n (inclusive). Set v i = 1 (arbitrarily), and 
use Equation[6]to obtain v 2 , ■ ■ . , Vd+i by setting 

(n - 2d)v i+1 - m 

v i+ 2 = ; (7) 

n — % 

for i = 0, 1, . . . , d — 1. Let p(x) be the unique polynomial of degree at most d 
with values p(i) = Vi for i £ {1, 2, . . . , d + 1}. Then the equation 

xp(x) + {n — x)p(x + 2) = Xp(x + 1) (8) 

holds, with A = n — 2d, for at least the d distinct values 0, 1, 2, . . . , d — 1 of x. 
Now both sides of Equation [8] are polynomials of degree at most d. If deg(p) < d, 
then both sides have degree less than d, so Equation [8] is an identity for p. If 
deg(p) = d, then both sides of Equation [8] not only have degree d, but also have 
the same leading coefficient of [n — 2d)cd, where q is the leading coefficient of p; 
subtracting this term from both sides, we find two polynomials of degree at most 
d — 1 which agree at d distinct values of x, so again they are identically equal. 
It follows that [p(l) p(2) ■ ■ ■ p(n + l)] 1 is an eigenvector of S n with eigenvalue 
n-2d. 

Since S n is an n + 1 by n + 1 matrix, it follows that the characteristic polyno- 
mial of S n is 

n 

Xn(x) = l[(x - (n - 2d)) = x £ - J] (x 2 - f), (9) 

d=0 l<]<n 

j—n (mod 2) 

lF The author wishes to thank R. Brualdi for informing him of this earlier work. 
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where e = if n = 1 (mod 2) and e = 1 if n = (mod 2). 

From Equation [9j we see that the coefficient of x n+1 ~ 2fc in Xn{ x ) ls 

d k , n ={-if Yi xif < io > 

H l t ] Y jeJ 

Now Equation Q] follows by comparing Equations l3l and [TOl 
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